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Abstract 

The effective operator tr(G 3 ) is the only dimension-6 gluonic operator that cannot 
be related to four-quark operators. A peculiar property of this operator is that it does 
not contribute to two-jet production at hadron colliders, at the level of one operator 
insertion and leading-order in a s ; therefore we study its effects on three jet events. To 
calculate the helicity amplitudes induced by this operator we make extensive use of collinear 
factorization. We propose several ways of detecting the tr(G 3 ) signal, one of which exploits 
its non-trivial behavior under azimuthal rotations of two almost collinear jets. 
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1. Introduction 

The gluon self-coupling is perhaps the clearest manifestation of the non-abelian nature 
of QCD. Tests of QCD are therefore incomplete without quantitative tests of the gluonic 
sector. New physics which respects the symmetries of QCD can produce deviations from 
the QCD prediction for the gluon self-coupling at higher energies. The deviations can 
be parametrized by an effective Lagrangian including higher-dimension operators which 
describe the low-energy effects of the new physics. By calculating the effect of these 
operators on partonic scattering, one can set lower bounds on the characteristic energy 
scale, A, of the new physics in a manner independent of the details of the new physics. 

The use of higher dimension operators as probes for new physics was originally sug- 
gested for interactions of leptons and quarks [1]. In the quark case, the dimension-6 op- 
erators are four-quark contact operators and they affect 2^2 quark scattering (qq —> qq 
etc.) at leading-order. Measurements of the dijet invariant mass distribution at 1.8 TeV 
by the CDF collaboration have recently led to a bound of 1.3 TeV on the scale A associated 
with such operators [2], where A is defined by the conventions of ref. [1]. 

This approach was first applied to testing purely gluonic interactions by Simmons in 
ref. [3]. It turns out, for reasons we will discuss in the following, that placing bounds on the 
scale of deviations from QCD in the gluon sector is more difficult than in the quark sector. 
It is this difficulty which the present work addresses. Though most previous work has 
concentrated on the two-jet cross-section, we will argue here that the gluonic dimension-6 
operator tr(G 3 ) can be bounded most efficiently by studying three jet events at hadron 
colliders. 

To test the gluon sector of QCD quantitatively, an effective Lagrangian that preserves 
the main features of QCD is desirable. At the level of dimension-6 operators the unique 
effective Lagrangian that can be constructed from gluon fields alone, and that is SU(3) 
gauge-invariant and CP-even is [3] 

Ces = £qcd + cj^fabcGTC^G^ + C'^D^D.G^ , (1) 

where G^ v = —{i/g s )[D^^ D v ] is the gluon field strength, D M = d^ + iggA^ is the covariant 
derivative, A and A' are the characteristic energy scales of the new physics and C and C 
are numerical coefficients. In accord with the conventions of ref. [1], we take C = C = An; 
then this equation defines A and A'. 

The new physics could be associated with heavy colored particles which couple to 
gluons and affect their interactions through loops, in which case the scales A, A' would 
be proportional to the mass of the heavy particle, with a proportionality factor that de- 
pends on the details of the new physics, such as the spin and color of the heavy particle. 
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Alternatively, the new physics could be gluon compositeness, in which case A, A' would 
be proportional to the scale of the interaction between the gluon constituents. In this 
work we do not concern ourselves with the possible source of new physics, but rather dis- 
cuss the implications of the effective Lagrangian (1) as a model-independent test of gluon 
self-interactions. 

Using the equation of motion for the gluon field A M , the second operator in (1) can 
be written as a four-quark contact operator [3]. Thus, it cannot provide a clean test of 
the gluon sector independently of quark effects. Also, since this operator is equivalent to a 
four-quark operator, it is fairly easy to obtain a good bound on the scale A associated with 
it just as in the case of four-quark operators; Simmons and Cho have recently estimated 
this bound as 2.03 TeV [4]. 

In the following we will therefore neglect the second operator in (1) (i.e., we set 
A' = oo), as well as four-quark operators, and focus on the first operator in (1), to which 
we refer as tr(G 3 ). 

Unfortunately, unlike the four-quark contact operators which interfere with pure QCD 
at leading-order in ot s and in 1/A 2 , the leading-order contribution of tr(G 3 ) to massless- 
parton scattering vanishes: the four-parton tree amplitudes containing one insertion of 
this operator have a helicity structure which is orthogonal to that of four-parton QCD 
tree amplitudes when the four partons are massless [3]. A leading order contribution 
from tr(G 3 ) does appear in four jet decays of the Z-boson, as discussed by Duff and 
Zeppenfeld [5], but because of their small center-of-mass energy these processes are only 
sensitive to scales A around 100 GeV, which is low compared to the typical energy scales 
of scattering processes in hadron colliders. Therefore, in order to obtain good bounds on 
the scale of tr(G 3 ) it is necessary to compute its effect on processes in hadron colliders at 
the next order, either in 1/A 2 or in a s . 

The former option was pursued by Simmons [3], and later by Dreiner et al. [6]. Sim- 
mons studied the effect of the gluonic operator on 2 — > 2 parton scattering (gg — > gg, 
QQ ~ *• 99 1 QQ ~ *• QQ an d processes related by crossing symmetry) at order 1/A 4 . Such 
contributions come either from squaring amplitudes with one insertion of tr(G 3 ), or from 
interfering QCD amplitudes with amplitudes containing two insertions of tr(G 3 ). The 
main feature of this correction as seen in dijet production, or alternatively, in inclusive 
jet production in hadron colliders, is an excess of events in the high p T region of the 
jet transverse-momentum (p T ) distribution. This signal roughly grows as (s/A 2 ) where 
s is the parton center-of-mass energy. It is therefore highly suppressed at low energies, 
and becomes appreciable only at energies of the order of A, where unitarity is violated, 
and where the parametrization in terms of the effective operator tr(G 3 ) can no longer be 
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trusted [6]. Including a form-factor to unitarize the cross-section reduces the tr(G 3 ) signal 
at high energies dramatically, leading to deviations from QCD that are on the order of the 
theoretical and experimental uncertainties [6]. 

An additional drawback of this approach is that at the order 1/A 4 contributions from 
dimension-8 operators must be included as well. Since there are many such dimension-8 
operators [3], it is very hard to treat the problem in full generality. 

The difficulties encountered here are related to the fact that the contribution studied 
is 0(1/ A 4 ). It seems desirable therefore to investigate the leading-order contribution in 
1/A 2 , while going to higher order in a s . In this case there are several options to consider. 
The first is to study the effect of the gluonic operator on inclusive jet production in hadron 
colliders at order 1/A 2 and next-to-leading order in a s . Three types of contributions appear 
at this order. The first is the interference of 2 — > 2 tree amplitudes containing one insertion 
of the gluonic operator, with 2^2 QCD loop amplitudes. This contribution is easy to 
obtain since the relevant QCD loop helicity amplitudes have been calculated using both 
string techniques [7] and Feynman diagrams [8] , and the helicity amplitudes containing one 
insertion of the gluonic operator are readily calculable with the use of the helicity basis. 
The second contribution involves the interference of QCD 2^2 tree amplitudes with 
loop amplitudes containing one insertion of the gluonic operator. The latter amplitudes 
are hard to compute. Finally, the third contribution arises from the interference of QCD 
2^3 tree amplitudes (gg — > ggg, gg — > qqg, qq — > qqg and processes related by crossing 
symmetry) with 2 — > 3 tree amplitudes containing one insertion of the gluonic operator. 

In addition to the technical difficulty of calculating the second contribution, this ap- 
proach suffers from a more serious problem. A rough estimate of the order of magnitude 
of (i) the QCD cross-section, (ii) the correction calculated by Simmons (0(1/A 4 )), and 
(iii) the (3(1/A 2 ) correction suggested above, reveals that the (9(1/A 2 ) correction only 
dominates over the (9(1/A 4 ) correction at low energies, where they are both insignificant 
compared to QCD. This is largely due to a factor of a s /(4-7r) suppressing the loop ampli- 
tudes. Instead, in this work we use the third contribution mentioned above, namely the 
interference of QCD 2^3 tree amplitudes with 2-^3 tree amplitudes containing one 
insertion of the gluonic operator, in order to probe the gluonic operator through exclusive 
three jet production in hadron colliders. In three-jet production, the tr(G 3 ) correction is a 
leading order correction, both in 1/A 2 and in a s . 

The use of three jet events in hadron colliders to test QCD has been discussed in the 
past [9,10]. Here we propose three different ways for studying the tr(G 3 ) signal in three jet 
events. The first involves the region in which the three jets are well separated, the second 
involves the region in which two of the jets are almost collinear, and the third involves the 
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cross-over between these two regions. The first approach is probably the easiest one from 
the point of view of the experiment but it leads to smaller tr(G 3 ) signals. 

The second approach is perhaps the most interesting from a theoretical point of view. 
In the region in which two of the partons are almost collinear, we use a special angular 
dependence of the tr(G 3 ) signal to isolate it from QCD. We discuss this approach in detail 
in section 3 and section 4 but we briefly outline the idea here. When two of the partons 
are almost collinear the three-jet event looks almost like a two-jet event, where one of the 
partons has a momentum p which is the sum of the two collinear momenta. The QCD cross- 
section and the tr(G 3 ) correction to it have different properties under azimuthal rotations 
of the collinear momenta around p, with p held fixed. An azimuthal dependence of the 
cross-section in the collinear region requires a linear polarization of the "effective gluon" 
which replaces the two collinear gluons in the resulting four-parton cross-section. However, 
a consequence of the supersymmetric Ward identities (SWI) [11,12,13] is that two of the 
three independent four-gluon helicity amplitudes, and one of the two independent two- 
quark two-gluon amplitudes vanish for tree-level QCD. It follows that for fixed helicities 
of the five external partons there is only one possible effective gluon helicity; it cannot be 
linearly polarized, and so the tree-level QCD cross-section has no dependence on azimuthal 
rotations of two collinear partons. On the other hand, after adding the tr(G 3 )-induced 
amplitude, both helicities contribute, giving rise to a linear polarization of the effective 
gluon and a characteristic azimuthal dependence of the cross-section. Precisely the same 
orthogonality of the helicity structures of the QCD and the tr(G 3 )-induced four-parton 
amplitudes, which prevents the appearance of the tr(G 3 ) signal at leading-order, also causes 
these different behaviors in the collinear region, thus providing a means for separating the 
tr(G 3 ) signal from the QCD background. In measuring the azimuthal dependence of the 
collinear three-jet cross-section, one simultaneously probes for the existence of the tr(G 3 ) 
operator, and tests the helicity structure of QCD and thereby the supersymmetric Ward 
identities. 

As is well known, the calculation of QCD amplitudes is not a trivial task. This 
task becomes even more daunting with the inclusion of the gluonic dimension-6 operator, 
which leads to cumbersome vertices for three-gluon, four-gluon and five-gluon couplings, 
all of which would be needed here. It is therefore crucial to use methods which simplify the 
calculation. This work illustrates the power of several such tools, namely, color ordering [14] 
the helicity basis [15, 16], and collinear factorization [17,16] for this extension of QCD. Thus 
for example, the orthogonality of QCD amplitudes and amplitudes induced by the gluonic 
operator at lowest order is immediately derived and becomes very transparent using the 
helicity basis. Calculating the new five-gluon vertex becomes unnecessary since we can infer 
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the five-parton amplitudes from the four-parton amplitudes using collinear factorization. 
For some of the calculations we can extract the contribution of the effective operator by 
thinking of it as induced by a heavy colored particle in the loop, taking advantage of 
existing string-based QCD loop calculations [7,18]. 

This paper is organized as follows: in section 2 we explain our notation and discuss the 
calculation of the new amplitudes. We give analytic expressions for the new four-parton 
and five-parton amplitudes, and for the corrections to the cross-section at the parton level. 
We discuss the behavior of these corrections in the soft and collinear regions of phase space 
in section 3. We then go on to describe the hadron-level calculation, and suggest different 
ways of looking for the tr(G 3 ) signal in three jet events in section 4. We conclude with 
a summary of the results in section 5. The appendix reviews properties of color-ordered 
amplitudes and collinear factorization, and shows that tr(G 3 ) does not induce any collinear 
singularities on top of the QCD singularities. 

2. Calculating the new amplitudes and interference terms 

The operator tr(G 3 ) induces new three-, four-, five- and six-gluon vertices. The three- 
gluon vertex has the same power of the strong coupling, g s , as the QCD vertex, and the 
other new vertices acquire an additional factor g s with each additional gluon. All the 
new vertices have 1/A 2 multiplying them. To compute the effect of tr(G 3 ) on three-jet 
production in hadron-colliders, one needs to evaluate five-parton tree amplitudes with one 
insertion of the new operator, i.e., one new three-, four- or five-gluon vertex. We denote 
these new amplitudes by the super-script (A). 

We now briefly describe our notation and the methods used to organize the calculation: 
color ordering, the helicity basis and collinear factorization. In this we follow closely the 
review by Mangano and Parke [16]. 

Tree amplitudes of SU(N) gauge theories can be organized in a simple color basis [16]. 
Each amplitude is written as a sum of color factors multiplied by color sub-amplitudes 
which correspond to specific color-orderings. Since all the color dependence is factored 
out of the sub-amplitudes, and since the color basis is orthogonal at leading order in the 
number of colors, each sub-amplitude is separately gauge invariant and can be calculated 
using any gauge choice. The n-gluon amplitude A n can be written as, 

An = Yl tr(T ai T a2 ---T a ") A n (l,2,.,n) . (2) 

{12. .n}' 

Here A n (l, 2, .., n) = A n (pi, \i;p2, A2; • • • ;p n , A n ) and pi, Xi are the momentum and helic- 
ity of the i-th gluon. The sum in (2) is over non-cyclic permutations of {1,2, ..,n}. We 
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discuss some useful properties of the color sub-amplitudes in the appendix. 
Amplitudes with n gluons and one quark pair can be written as, 



A n = ^ (T a *T a *...T a ") id A n (q ] l,2,..,n ] $. ( 3 ) 
{12. .n} 

Here A n (q; 1, 2, .., n; q) = A n (q x ;p 1 , Ai;p 2 , A 2 ; • • • ;p n , A n ; q X ), q, A are the quark momen- 
tum and helicity, q, A are the anti-quark momentum and helicity, Pi, Aj are the momentum 
and helicity of the i-th gluon, and the sum is over all permutations of the n gluons. 

We will also need amplitudes with two quark pairs and one gluon. These we write as, 

As = A 5 (q 1 ,q 1 ;q 2 ,q 2 ;k) -j^S^T^ + A 5 (q u q 2 ; qz, qi\ k) S ilia T^^ + 

1 " " ' (4) 

+ A 5 (q 2 ,q2;qi,qi;k) —S^T?^ + A 5 (q 2 , qi; q u q 2 ; k) S^T^ , 

where qi, q~i are the momenta of the i-th quark and anti-quark, and k is the gluon momen- 
tum. 

Calculating the new amplitudes is greatly simplified by using the helicity basis [15,16]. 
Each gluon polarization is represented in terms of the gluon momentum and another light- 
like momentum, the reference momentum. Choosing the reference momentum as an ap- 
propriate combination of the external momenta (this choice can be made separately for 
different color sub-amplitudes) results in many cancellations in the calculation. With this 
representation of the polarization vectors, helicity amplitudes of massless external partons 
are naturally expressed in terms of spinor products, which are nothing but square roots of 
the Lorentz invariants up to phases. The spinor products are denoted as: 

[ij] = (i+\3~) = u + (ki)u-(kj) , 
{ij) = (i-\j+) = u-ikiju+ikj) , 

where )) is a positive (negative) helicity spinor with the momentum of parton i. 

We give an explicit representation of the spinor-products in the appendix. 

One of the useful properties of color-ordered helicity amplitudes is their factorization 
on collinear poles. When the color adjacent partons i, i + 1, with momenta ki, ki+i 
become collinear, the n-parton color sub-amplitude factorizes on the (n — l)-parton color 
sub-amplitudes where the two collinear partons are replaced by one parton with momentum 
p = ki + ki + i and with positive or negative helicity. We sometimes refer to this parton as 
the effective parton. Except for a subtlety in the terms sub-leading in the number of colors 
in the four-quark one-gluon amplitudes (discussed below), each n-parton color structure 
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factorizes on a unique (n — l)-parton color structure. In the collinear limit ki = zp, 
ki + i = (1 — z)p and [17,16], 



A n (k u ...,ki, k i+1 , ...,k n ) ^ t+1 Split_ A (z, % + 1) A n _i(fci, . . . , (ki+k i+1 ) x , ...,k n ). 

(6) 

The splitting function Split_ A (i, j) depends on the helicity A of the effective parton, on 
the helicities of the two collinear partons and on the momentum fraction z, and is singular 
as s^i+i = (ki + fci+i) 2 — > 0. The splitting functions are analogues of the Altarelli- 
Parisi coefficients [19] that are associated with color-ordered helicity amplitudes instead 
of cross-sections. Specifically, the absolute value of the splitting-function, squared and 
color-averaged, is equal to the appropriate polarized Altarelli-Parisi coefficient. We list 
the relevant splitting functions in the appendix. 

Amplitudes with insertions of tr(G 3 ) satisfy (6) as well. When the two collinear 
gluons z, j are attached to the new three-gluon vertex, one can show (see appendix) that 
the expression for this vertex goes to zero at least as fast as s^-, cancelling the 1/s^ 
propagator pole. This is related to the fact that the new vertex has extra powers of 
momenta in the numerator compared to the QCD vertex, to compensate for 1/A 2 . Thus 
the only splitting functions which occur for the five-parton new amplitudes are the standard 
splitting functions arising from the QCD vertices and so, 

A^\k u ...,ki, k i+1 , ...,k n ) J2 SpHt_ A (i, i + 1) • • • > (h+k l+1 )\ . . . , k n ) . 

A=± 

(7) 

Collinear factorization provides a check of n-parton amplitudes when the (n — 1)- 
parton amplitudes are known. Moreover, it is often possible to infer the former from the 
latter [20]. In this work we calculate the relevant tr(G r3 )-corrected four-parton amplitudes 
and then guess the form of the five-parton amplitudes by requiring that they have the 
correct behavior in all collinear limits. There could be, however, additional terms in the 
amplitudes which are finite in all the collinear limits. We have not been able to find 
any such terms that also have the correct dimension and for the five-gluon amplitudes also 
satisfy the U(l) decoupling equation (see the appendix). Still, to verify the results obtained 
from collinear factorization, we calculated at least one helicity amplitude for each case (5- 
gluon, 2-quarks 3-gluons, and 4-quarks 1-gluon) using an alternative method. We find 
agreement with the results obtained from collinear factorization; there are no additional 
finite terms which are missed when collinear factorization is used to infer the five-parton 
tr(G 3 ) amplitudes. 
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We begin by listing the relevant four-parton amplitudes. These are obtained from 
Feynman diagrams using the helicity basis, and can be written as one-term expressions. 
For the four helicity amplitudes in gg — > gg we have, with g s = 1, 

4 A) (l+,2+,3+,4+) = * 2gl2g23g13 



A 2 (12) (2 3) (3 4) (41) ' 

2ro Al2r, r, n 2 



A (A )( + + +) _ 3z -[23]-[34]-[42]- 
Aa [L ,Z ,<5 ,4 J - A2 [12] [2 3] [34] [41] ' ( 8 ) 



4 A) (l-,2-,3+,4+) = 0, 
4 A) (l-,2+,3",4+) = 0. 

The corresponding gg — > gg QCD tree amplitudes are: 

A 4 (l+,2+,3+,4+) = 0, 
44(1-2+3+4+) = 0, 

(34) 4 



(9) 



A 4 (l+,2+,3-,4-) = A 4 (l+,2-,3+,4-) 



(12) (2 3) (34) (41) 



Comparing (8) and (9) it is easy to see that the interference vanishes, and there is no 
tr(G 3 ) induced correction at leading order. Similarly, for qg — > qg, 



4< A >f*+-1 + 2+-a-) 3z [12]-[gl] 2 [Q2] [2q] 

4 {q ' ' ' q } ~ A 2 [gg][gl][12][2g] ' 

A^(a~-1 + 2+-o+) - ^ -[12] 2 [g2] 2 [ 9 l] [lg] 

4 ( " ' ' } " A 2 [g 3 ][gl][12][2g] ' ( 10 ) 

4 A) ( g +;l+,2-;r) =0, 

4 A V;1+,2+;q+) = 0. 

The g(7 — > QCD tree amplitudes are: 

A 4 ((z+;l+,2+;n = 0, 

A 4 (q-;1+,2+;q+) = 0, 

AJa+-l 2-a~) - i (11) 

Mm,*,?) = * {qifiqi) 



(qq) (ql) (12) (2q) 



where i is the negative helicity gluon and again there is no leading-order interference. 

An alternative way to obtain the amplitudes in (8) and (10), is to think of the operator 
tr(G 3 ) as being induced by a heavy colored scalar of mass M s on the order of A, circulating 
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in a loop, and to calculate the 1/M S 2 term in the resulting loop amplitude. This is not as 
cumbersome as it may seem, since all the relevant loop amplitudes have been calculated 
by Bern and Kosower [7] using string techniques for the case of a massless particle in the 
loop. The Feynman parameter polynomials for a massive scalar are identical to those for a 
massless scalar, since the derivative interactions are identical. Only the scalar denominator 
of the Feynman integral is changed. In the limit of large M s , the only (3(1/M 2 ) contribu- 
tions come from the triangle diagrams, i.e., diagrams with three legs attached to the loop. 
The resulting Feynman parameter integrations are simply integrations over polynomials. 
Subtleties of ultraviolet and infrared divergences and coupling-constant shifts, do not oc- 
cur in these loop amplitudes because the corresponding QCD tree amplitudes vanish. By 
calculating A^\q + ; 3 + , 4+; q~) in both ways, we extract the ratio between M s and A, 
A 2 = 7207rM 2 /a s . We then use this ratio to check that the two methods agree for the 
remaining amplitudes. 

The zeroes in equations (9) and (11) are manifestations of the supersymmetric Ward 
identities (SWI) [11,12] (unlike the zeroes in equations (8) and (10), whose origin is ob- 
scure). According to these identities, the only non- vanishing four-parton amplitudes are 
the ones with two positive and two negative parton helicities. For QCD, which is not 
supersymmetric, the SWI are satisfied at tree-level only [13]. Because of the SWI, as we 
will see in the following, the five-parton tree-level QCD cross-section is invariant under 
azimuthal rotations of two collinear partons which leave the sum of the collinear momenta 
fixed. 

We now turn our attention to the five-parton amplitudes. Since we are interested 
in the interference of tr(G 3 )-induced amplitudes with QCD tree amplitudes, we have to 
calculate only the new amplitudes corresponding to helicity choices for which the QCD 

tree amplitudes are non-zero. These are (+ + H ) and (+ H 1 — ) for the five-gluon 

amplitudes, (±; + H — ; =f), (±; H h; =f) and (±; — h +; =f) for the two-quark-three-gluon 

amplitudes, and all possible helicity structures for the four-quark-one-gluon amplitudes. 

To illustrate the use of collinear factorization in deriving these amplitudes, consider 
the 2-quark 3-gluon amplitude A^ A \q + ; 1 + , 2 + , 3~; q~). In the limit when q and k\ become 
collinear, this amplitude should factorize on the 4-parton amplitude A^\p + ; 2 + , 3~; q~) 
with p = q + k\. Since the latter amplitude is zero, Al A \q + ; 1 + , 2 + , 3~; q~) cannot be 
singular in this limit. Similarly, it remains finite as gluons 1 and 2 become collinear 
because Split + (+, +) is zero (see the appendix). However, in the limit where gluons 2 and 
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3 are collinear with k 2 = zp and ks = (1 — z)p, 

4 A) (q+;1+,2+,3-;q-) ~ Split. (2+ 3") x A[ A) (q+ ] l+,p+; q~) 

- i/2 3 .1 [ip][i?][pg] 



(l-z) 1/2 [2 3] A 2 [9 9] 
This is satisfied by the expression 

A^(a + -1 + 2+ 3~- a') - 3i- MMMM 
A 5 (q ,1 ,2 ,3 ,q ) - 3i ^ { - q] [2 3] [s - } , 

which also has the correct behavior when g and £3, and q and g become collinear, and 
which agrees with the result obtained from a Feynman diagram calculation. 

In tree amplitudes with zero or one quark pair, the n-parton color sub-amplitude 
always factorizes on a unique (n — l)-parton sub-amplitude. The situation is slightly 
different for amplitudes with two quark pairs, where the color ordering is not unique 
at sub leading order in the number of colors. Thus for example as q\ and q 2 become 
collinear in the 0(1/N C ) sub-amplitude A( A \qi, q~i; q 2l q 2 ; k) there are contributions from 
both orderings of the two gluons (the original gluon and the one replacing the pair qi, q~i) 
attached to the remaining quark line. Thus one finds, 

A {A \qi,qi,q 2 ,q 2 ;k) ~ Split + (gi, qi) (A(q 2 ; k,p; q 2 ) + A(q 2 ; p, k; q 2 )) . 

Our results for the tr(G 3 ) -induced five-parton amplitudes are as follows. 
For gg -> ggg, 

A {A) (i+ 7+ k+ r m-) ~ - ^ ij ^ k ^ ki ^ (12) 
A 5 [i ,j ,fc ,1 ,m ) - A2[1 2 ] [ 2 3] [3 4] [4 5] [51] ' [ ] 

where l~ and m~ are not necessarily adjacent. This expression was inferred from collinear 
factorization. As a check, we modified the string-based calculations of five-gluon loop 
amplitudes with a massless scalar in the loop [18] to obtain the 1/Mf contribution from a 
scalar of mass M s , which agreed with equation (12). 

The gg — > ggg QCD tree amplitude which interferes with (12) is: 

A 5 (i+, j\ r,m-) = i (i2) (23) ( | 3 m 4 ) )(45)(51) > ( 13 ) 

where, as in (12), / and m are not necessarily adjacent. 
For the two-quark three-gluon amplitudes we find, 

5 {q ' ' ' 3 '" } " A 2 [^][al][12][2 3][ 3(? ] • ( 14 ) 
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In this formula i and j are the two positive helicity gluons, and i is closer to q than j is in 
the color ordering. The formula for reversed fermion helicities is 

5 {Q ' 1 ' 2 ' 3 ' 9 } - AM^][ ? 1][12][2 3][3,1 ' (15) 

where now z is the gluon closer to These amplitudes were constructed to have the 
correct collinear limits. As a check we calculated A^ A \q~; 1 + , 2 + , 3~; q + ) from Feynman 
diagrams. 

Notice that the amplitude with reversed quark pair helicities can be obtained by 
exchanging the quark and anti-quark and reversing the ordering of gluons along the quark 
line; changing the helicities of the quark-pair amounts to a charge conjugation operation 
which reverses the color ordering associated with the quark line. Thus equation (15) can 
be obtained from equation (14). 

The required QCD amplitudes are: 



A,(a+-1 2 3-0-1 - i ZMM! 
5( " ' 1,2,3,9 } " \qq)(ql)(12)(23)(3q) 

Mq-,l,2,3;t) = i {qi){qif 



(16) 



(qq) (ql) (12) (2 3) (Sq) 



where i is the negative helicity gluon. 

Finally, for the four-quark-one-gluon amplitudes we have, 

M A )(a+ a--a+ a--k+) - 3z [gj ggj [gj k ^ [gg ^ 

A 2 [?l?l][?2?2] (17) 

A (A) (qi,q i; q 2 ,q 2 \k) = . 

These amplitudes were inferred from collinear limits, and verified by a Feynman-diagram 
calculation. Again, amplitudes with reversed quark pair helicities can be obtained by 
exchanging the quark and anti-quark, and this can be done for each quark-pair separately. 
The corresponding QCD amplitudes are: 

A(qt^;qt^k + ) = F {<l2<i:) 



(q2k){ qi k) ' 

A(qt,q^,q+,q^,k+) = F 



(q 2 k) (q 2 k) 
where 

F = i ^ 



(Qi Qi) (92 Q2) 
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Amplitudes with reversed quark-pair helicities can be obtained from (18) by exchanging 
the appropriate quark and anti-quark in F. 

Notice that the expressions in (17), (18) were derived for different quark pair flavors. 
The amplitudes with same flavor quarks can be obtained from them by antisymmetrizing 
over identical quarks. 

For all the different parton processes, one can obtain the amplitudes with all helicities 
reversed by parity, which exchanges angle-brackets ( ) and square-brackets [ ] in (12) - (17), 
and multiplying by an additional minus sign for amplitudes with an odd number of gluons. 

We now need to interfere these amplitudes with the appropriate QCD amplitudes, 
and to sum over colors and helicities. Since the initial and final states vary from process 
to process, we perform the initial state averaging later. 

For gg — > ggg we get, summing over final states colors and helicities, 



UU 99999 



helicities colors 

= -6^ 6 ^iV c 3 (iV c 2 - 1) x 



)- 



x 



([12] [23] [34] [45] [5 I])" 2 £ ([ij] [j k] [ki] [I mf) + c.c. 

{2,3,4,5} Z,m6{l,2..,5}^ m 



(19) 



where the first summation is over all permutations of (2, 3,4, 5). 

For scattering processes involving four quarks and one gluon, the interference term 
summed over colors and helicities is, 



™ ^A 2 C[ ° (q iqi ) 2 (q 2q2 f Ul *> <<fc *> <&*><«!*> 



x 



x (Qi Q2) 2 (qi ^2) (qi k) (q 2 k) + (qi q 2 ) 2 (<?i 92) (<?i k) (q 2 k) + (20) 



+ {Qi Q2} 2 (qi 92) (qi k) (q 2 k) + (qt q 2 f (q x q 2 ) (q x k) (q 2 k) + c.c. 



Again, equation (20) holds for different quark pair flavors, but the modification for same 
quark flavors is straightforward. 

The corresponding expressions for scattering processes involving two quarks and three 
gluons are more complicated, so we only give here the analytic expression for the color 
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sum, without summing over helicities. (We perform the latter summation numerically.) 



&» = %^ E f^ 2 ) 2 E ^'(1,2,3)^(1,2,3) + 



fo " 

' helicities 



{1,2,3} 



+ iV c 2 A (A) (1,2,3) (-2A*(1,2,3) - A*(2,l,3) - A*(l,3,2) + 



{1,2,3} 



(21) 



+ A*(3,2,1)) + ( ^ (A) (1,2,3))( ^(1,2,3)) + c.c. 



{1,2,3} {1,2,3} 

Here all summations are over the six permutations of (1, 2, 3), and A(l, 2, 3) is shorthand 
for A(g;l,2,3;g). 

3. Collinear and soft behavior of the tr(G 3 ) correction to the cross-section 

As is well known, the leading order QCD cross-section for three jet production is 
singular in the soft and collinear regions of phase-space. This property is exhibited as 
strong peaking of the QCD distributions near values of the three-jet event variables which 
correspond to a collinear or to a soft region. In order to search for the tr(G 3 ) signal in three 
jet events it is important to understand its behavior in these regions. As mentioned in the 
introduction, the special angular dependence of the tr(G 3 ) signal in collinear regions can 
in principle be used to separate it from QCD. In this section we will explain this angular 
dependence in more detail. First, however, we discuss the behavior of the tr(G 3 ) signal in 
soft regions. 

Inspecting the expressions for the interference of QCD amplitudes and tr(G 3 )-induced 
amplitudes reveals that the different parton energies appear with the same power in the 
numerator and the denominator in each term, so that the interference remains finite as 
any of the energies approaches zero. Thus the tr(G 3 ) correction to the three jet cross- 
section remains finite in soft regions of phase-space. Appropriate cuts should therefore be 
imposed to avoid these regions in searches for tr(G 3 ) since the singular QCD background 
will swamp the signal there. 

As for collinear regions, we know that the new amplitudes have collinear singularities 
since these were used to construct them. For both the QCD and the new amplitudes, the 
singularity associated with partons i,j becoming collinear is of the form (ij) 1 or 1 , 
so that its magnitude is Sy _ 3. Therefore it results in an Sjj -1 singularity in the QCD 
cross-section. 

In the following collinear analysis, we hold fixed the helicities of the five external 
partons, and examine the possible helicities of the effective parton with momentum ki + kj. 



14 



If only one of Split + and Split_(z,j) is non-zero, the resulting singularity in the 
tr(G 3 )-corrected cross-section is only Sjj - ^, because the helicity structures of the four- 
parton QCD amplitudes and the new amplitudes are orthogonal, so that only one of the 
two amplitudes in the interference is singular. However, for gg — > g and qq — > g both 
Split + (z, j) and Split_ (i, j) can be non-zero, in which case the tr(G 3 ) correction to the 
cross-section has a single pole in Sy. Then, the QCD cross-section and the new cross- 
section are equally singular in the limit when two partons become collinear. Still, the 
form of the singularity is different. The singular term in the tree level five-parton QCD 
cross-section is always of the form |Split ± (z, j)\ because there is only one choice of the 
effective parton helicity that yields a non-zero four-parton amplitude, due to the SWI 
zeroes in equations (9) and (11). In contrast, in the interference of a QCD amplitude 
with a tr(G 3 ) amplitude, the choice of the effective parton helicity is still unique in each 
amplitude, but it must be of the opposite helicity for the QCD amplitude as for the tr(G 3 ) 
amplitude, due again to four-parton orthogonality. Therefore, the singularity is always of 
the form Split ± (i, j)*Split =p (i, j). This last factor carries the phase 2(p, where (p is the angle 
associated with azimuthal rotations of the two collinear momenta around the direction of 
their sum, p, with p held fixed. Notice that the phase has to be 2ip rather than ip or it 
would matter which of the two partons i, j we use to determine <p (exchanging % and j 
amounts to exchanging (p and it — (p) . This is the only (p dependence of the cross-section 
in the collinear region, as the rest of the expression is just the product of the QCD and 
tr(G 3 ) four-parton amplitudes and they only depend on the four-parton parameters. 

As we mentioned above, the trivial p- dependence of the QCD cross-section follows 
from the supersymmetric Ward identities. The QCD cross-section has no p> dependence 
because of the zeroes in equations (9) and (11), and these, in turn, are a consequence of 
the SWI. 

The ip dependence of the cross-section is related to the polarization state of the effec- 
tive parton. Each QCD five-parton amplitude factorizes on just one four-parton amplitude 
with either a positive or a negative circularly-polarized parton. The resulting QCD cross- 
section can therefore have no tp dependence. In contrast, each tr(G 3 )-corrected five-parton 
amplitude factorizes on the sum of two four-parton amplitudes with oppositely polar- 
ized gluons. The resulting cross-section contains an interference term corresponding to a 
Zmear%-polarized effective gluon, and this term gives rise to a non-trivial ip dependence, 
correlated with the polarization vector. 

Since the singular part of the tr(G 3 ) correction behaves as e 2lip it is washed out upon 
integrating p> between zero and 2rc. The tr(G 3 ) correction to the total cross-section is 
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therefore finite; it has neither soft nor collinear singularities f. 

The non-trivial <p dependence of the tr(G 3 ) signal can be used to separate it from 
tree level QCD. If we consider three jet events in a region in which two of the partons are 
almost-collinear, it is always possible to choose a frame in which the (p dependence of the 
tr(G 3 ) correction is cos2</?. If we then weight the events by cos2</? and integrate over (p, 
tree level QCD washes out while the tr(G 3 ) contribution remains non-zero. It is useful to 
define the cos 2<p expectation value, 



where M and M' are either the QCD or the tr(G 3 )-induced five-parton amplitudes (de- 
pending on whether one is considering the pure QCD signal or the tr(G 3 ) correction to the 
signal), the sum is over colors, helicities and the different relevant parton processes, and 
the integration is restricted to the almost-collinear region of phase space and sweeps the 
full range of ip, < ip < 2tt. This expectation value is a good probe of the operator tr(G7 3 ); 
it receives contributions from the tr(G 3 ) correction and no contribution from leading-order 
QCD. Furthermore, it provides a test of the SWI, or alternatively, of the helicity structure 
of tree-level QCD. 

Notice that four-quark contact operators do not contribute to (cos2</?) at leading- 
order. These operators only give rise to four-quark one-gluon amplitudes, which can only 
factorize on four-quark amplitudes containing the four-quark contact vertex, when a quark 
and a gluon become collinear. Since the helicity of the effective quark is determined by 
the helicity of the original quark, the same splitting function occurs in both these new 
amplitudes and the corresponding QCD amplitudes so that the interference has no ip 
dependence. Thus, these operators do not contribute to the cos 2ip expectation value. 

However, some background to the tr(G7 3 ) signal as measured through this expectation 
value will arise from higher order QCD effects. A similar azimuthal dependence to the 
one described above occurs in next-to-leading-order corrections to 2 — > 3 scattering and in 
tree-level 2^4 scattering in pure QCD. The first effect turns out to be very small. The 
collinear behavior of five parton loop amplitudes is similar to equation (6), except there are 
now two terms [21]. In the first term a "loop splitting-function" multiplies four-point tree 
amplitudes; here the previous tree-level arguments still forbid an azimuthal dependence. 

f This implies that the four-parton loop amplitudes with one insertion of the effective operator, which 
were mentioned in the introduction, are finite. There are neither ultraviolet nor infrared divergences 
(soft nor collinear) for them to cancel against. In addition, these loop diagrams have no imaginary 
parts, as can be seen by cutting them into two tree diagrams and using the orthogonal helicity 
structures of the QCD and tr(G 3 ) four-parton tree diagrams. So they are apt to be quite simple 
expressions. 




(22) 
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In the second term, the usual tree splitting-function multiplies a four-parton QCD loop 
amplitude; the only azimuthal dependence arises when the loop amplitude is one which 
vanishes at tree-level. The relative magnitude of this effect compared to the tr(G 3 ) signal 
is given by the relative magnitudes of the tr(G 3 ) correction to the four-parton amplitude 
and the relevant (infrared and ultraviolet finite) four-parton QCD one-loop amplitude. 
The relevant four-gluon QCD loop amplitude, for example, is [18] 

A i-ioo Pfl _ 2+ 3+ +) _ _J_ N L _ nA (2 4) [2 4] 3 

Aa U ' 2 ' 6 ' 4 ] ~ 48tt 2 ' c { N C J [12] (2 3) (3 4) [41] ' 

Including the coupling constants, it is smaller than the tr(G 3 ) amplitude (8) by a factor 
of the order of (A 2 /s)(a s /187r). For the values of A and s relevant for the present analysis 
(see section 4), this ratio is of the order of 10~ 3 — 10~ 2 . This estimate is essentially the 
same as the one mentioned in the introduction where it was remarked that the interference 
of a QCD loop amplitude with a tr(G 3 ) amplitude is small compared with the square of a 
tr(G 3 ) amplitude. A similar estimate holds for quark-gluon amplitudes. 

The second source of background, 2^4 scattering in pure QCD, is harder to estimate. 
In this case too there could be a non-trivial dependence on azimuthal rotations of two 
almost collinear partons, but to get a large correction to the three-jet cross-section one of 
the remaining two final state partons should be soft or collinear. However, in this limit the 
six-parton amplitude factorizes on the five-point tree amplitude so there is again a unique 
possibility for the helicity of the effective parton leading to the conclusion that 2^4 QCD 
contributions to (cos 2(p) are probably small. 

4. The tr(G 3 ) signal in three jet events 

We can now proceed to calculate the effect of tr(G 3 ) on different distributions in 
three jet variables. We specifically consider here three jet production at the Tevatron as an 
example. Since we are only interested in a qualitative estimate of the tr(G 3 ) signal, we take 
the jet momenta to coincide with the parton momenta, without using any hadronization 
algorithm. We still have to fold the partonic cross-section with the proton, or anti-proton, 
structure functions: 



da 



J dx 1 dx 2 F 1 (x 1 ,Q 2 )F 2 (x2,Q 2 )da i , (23) 



where i labels the different partonic processes, and we choose the scale Q as Q 2 = {x\X2s). 
We use the DFLM [22] structure functions for the calculation. Here dai stands for either 
the QCD cross-section or the tr(G 3 ) correction. We evaluate the two simultaneously, 
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which allows for some checks on the calculation by comparing our results for pure QCD 
to existing simulations of the QCD distributions (we used the CDF collaboration's results 
as reported in [23]). Since there is no mass scale in the problem besides A (we neglect 
quark masses), both the QCD cross-section and the tr(G 3 ) correction have definite energy 
scaling properties; the QCD cross-section is energy-independent, and the correction scales 
as the energy squared. This can be used to factorize the calculation into two parts: an 
angular integration over final states, and an integration over the incident partons' x\ and 
x 2 . This leads to a significant reduction in computer time. We perform the first integral 
using the Monte-Carlo program SAGE [24]. 

We also take into account the renormalization of tr(G 3 ) between the scales A, where 
tr(G 3 ) appears with a coefficient of 4-7T (see equation (1)), and the scale Q [4]. The anoma- 
lous dimension of tr(G 3 ) is equal to \g 2 s /&n 2 with A = 7 + 2n//3 [25,26] so that the 
coefficient of tr(G 3 ) at the scale Q is, 



where b = —11/2 + n//3 and, 

a s (Q) 



a s (A) 



a s {M z 



x 

2b 



C (Q) = > ( 24 ) 



l- 6 £^log£ 



M z 

where we take a s {M z ) = 0.118. The resulting reduction in the coefficient of tr(G 3 ) for 
the relevant values of A and Q (see the following) is roughly ten percent. The operator 
tr(G 3 ) does not mix with any other dimension-6 operator through the renormalization 
group equation [25,26]. In this sense, studying it separately from the other operators is 
not inconsistent. 

We now discuss the geometry of three-jet events (see fig. 1), the cuts we impose and the 
different distributions we study. We use the notation of the CDF collaboration [23] . Viewed 
in the center-of-mass system, the three outgoing jets can be described by five independent 
variables; the total energy of the jets V§, and four scale-independent quantities, namely, 
the energy fractions of two of the jets, and two angles. Labeling the jets 3, 4 and 5, in 
order of decreasing energies, the energy fraction of jet % is defined as 

Xi = ^| , i = 3, 4, 5 , (25) 

vs 

where Ei is the energy of the i-th jet, and xs + X4 + x$ = 2. The angles we use are the 
angle between the fastest jet and the beam direction, 9, and the angle between the jets' 
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plane and the plane of the fastest jet and the beam, ip. An additional angle, the overall 
azimuthal angle around the beam direction, is required for a full description of the event, 
but the dependence on it is trivial for unpolarized beams. 

We describe here three different ways of probing tr(G 3 ). The first is to study distri- 
butions in various three jet variables in the region where the three jets are hard and well 
separated. The second is to study the cos 2ip expectation value of equation (22) obtained 
by weighting the relevant matrix elements by cos 2<p in the region where jets 4 and 5 are 
almost collinear as discussed in section 3. Finally, one can interpolate between these two 
regions and study distributions which reflect the qualitative differences between the QCD 
cross-section and the tr(G 3 ) correction in the cross-over region. 

To ensure that the three jets are well separated we impose the following set of cuts 
(set A) 



v 7 ! 


> 


250 GeV , 


x 3 


< 


0.8 , 


cos6> 


< 


0.8 , 


30° 


< 


%j) < 150° 



(26) 



We take the lower bound on the total jet energy, y/I, to be 250 GeV for the Tevatron 
following ref. [23]. Increasing this bound enhances the deviation from QCD, but decreases 
the cross-section, thus leading to poorer statistics. The energy fraction of the fastest jet, 
#3, varies between | and 1. The former value corresponds to a symmetric event, and the 
latter to either a soft jet (5) or to two collinear jets (4 and 5). It is therefore necessary 
to set some cut on the maximum value of x 3 to avoid soft or collinear jets. As this cut 
decreases, the sensitivity of the signal to the deviation from QCD improves, since it drives 
the events away from the soft region, where QCD is divergent whereas the deviation is 
finite. We take the x% cut to be 0.8. To avoid regions where jet 3 becomes collinear with 
the beam direction, we require that the absolute value of cos 9 be smaller than 0.8. Finally, 
to prevent jets 4 and 5 from approaching the beam direction we take to lie between 30° 
and 150°. 

The distributions in X3, ip and the three jet invariant mass a/S for the scale A = 1 TeV 
are plotted in Fig. 2. The errors bars shown are the estimated statistical errors assuming 
an integrated luminosity of 15pb _1 . The tr(G 3 ) correction to all three distributions is 
positive throughout this range. The deviations from QCD change sign because we plot 
normalized distributions following the CDF conventions. Even though the shape of the 
tr(G 3 ) correction is different from the shape of the QCD distributions for these variables, 
when added to the larger QCD cross-sections they still result in distributions that are 
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qualitatively similar to the QCD distributions, apart, roughly, from a vertical shift. Thus, 
the effect of tr(G 3 ) on the normalized distributions is quite small in this region, and 
although it is larger than the estimated statistical errors, it seems unlikely that studying 
this region by itself could conclusively determine the presence of the effective operator. 

We have also computed distributions in cos _1 ((x4 — x^)/xz), as suggested by Ellis, 
Karliner and Stirling [10], in (1 + cos 6>)/(l — cos (9), following ref. [9] and in X4 and 9. None 
of them improves the sensitivity to the signal. 

The situation becomes even worse if one includes a form-factor in the analysis. The 
operator tr(G 3 ) violates unitarity in 2 — > 2 parton scattering as discussed by Dreiner et 
al [6]. These authors replace 1/A 2 by 

To " 9 ( 27 ) 

A (1 + 47rs/(3A 2 )) 

where s is the parton center-of-mass energy, to restore unitarity in most partial waves in 
the 2^2 scattering f . Since the tr(G 3 ) signal in dijet production is (9(s 2 /A 4 ) it is greatly 
reduced by the form-factor. In the case at hand, the effect of the form-factor is smaller 
but it still reduces the already-small tr(G 3 ) signal of Fig. 2. 

Alternatively, the effective operator can be probed in the second way mentioned above, 
namely, by focusing on the collinear region and using the ip dependence of the correction 
induced by the effective operator to detect it. To get to the region where jets 4 and 5 are 
almost collinear and neither is soft we take xs > 0.95 and x§ > 0.3. We need some further 
cut to ensure that the two jets are still distinguishable. If we were to use a standard jet 

/ 2 2 

algorithm which involves R = y (At]) + (A^) , where rj is the jet's rapidity and 4> is the 
usual azimuthal angle of the jet, in order to acheive this, we would introduce a spurious 
(/^-dependence which would spoil the effect we describe here. Instead it would probably be 
best to impose a cut on the minimum transverse momentum of jets 4 and 5 with respect 
to the direction of jet 3 f (this direction is the direction of the boost to the rest frame of 
jets 4 and 5). However, since we factorize the calculation into independent integrations 
over xi, X2 and over the final states, we impose the stricter cut: 

• „ 5GeV 
£4 sin #34 > — — , (28) 

W S )miJ 2 

where 6*34 is the angle between jets 3 and 4. Combined with the cut on s, s > s m i n , this 
cut guarantees that the transverse momentum of jets 4 and 5 is greater than 5 GeV. 



f Here we modified the numerical factor to take into account the different scale definitions. The scale 

A used by Dreiner et al. is equal to our A/VI671-. 
t We thank J. D. Bjorken for suggesting this. 
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In addition, we take | cos6>| < 0.5 so that the 2-jet-like event is almost perpendicular 
to the beam direction and the tr(G 3 ) correction is more pronounced. However, we impose 
no cuts on the angle ip. This angle should be allowed to sweep its entire range so that the 
azimuthal angle <p of jets 4 and 5 can vary between zero and 2tv. 

We collect these cuts (set B) here 



v s 


> 250 GeV , 


x 3 


> 0.95 , 


x 5 


> 0.3 , 


X4 sin 6*34 


5GeV 


" (^) m ; n /2 


cos6> 


< 0.5 . 



(29) 



The resulting QCD and tr(G 3 ) -corrected distributions in x 3 and in the cosine of the 
angle between jets 4 and 5, cos #45, are plotted in Fig. 3. Again, the error bars shown are 
the estimated statistical errors assuming an integrated luminosity of 15pb _1 . In 3a and 
3b we plot the distributions obtained by weighting the matrix elements with cos 2(p, 

dotn f d 2 a da 

- 1 dip cos2</? — — — = (cos2</?) — , (30) 



d cos 6*45 J dtpd cos $45 d cos 9 45 

and, 

da^ f , „ d 2 a , „ , da 



. dip cos2(/9 - — - — = (cos2<^) - — , (31) 
dxs J dipdxs dxs 

where (cos2<^) of (30) is different than that of (31) as they are distributed in different 
variables. In this case the angle (p is the azimuthal angle of jet 4 around the sum of the 
momenta of jets 4 and 5, and it is related to the angle ip: <p = ip — ir/2. The resulting 
distributions are negative as they receive larger contributions from events that are co-planar 
with the beam. The magnitude of the QCD distributions decreases on approaching the 
collinear limit while the magnitude of the tr(G 3 ) correction grows. In 3c and 3d we show 
(cos2<^) as obtained from QCD with and without tr(G 3 ). Clearly, the signal is dominated 
by the tr(G s ) correction. Also shown are the tr(G 3 )-corrected distributions with the form- 
factor (27) included. The form- factor reduces the tr(G 3 ) correction uniformly over the 
entire x s and cos 6> 45 ranges. The reduction factor is equal to the scale-averaged form- 
factor, and is about one half, which is roughly the value of the form- factor near the cut 
= 250 GeV, from where most of the cross-section comes. 

Recall that the tr(G 3 ) correction to the cross-section has a non-trivial ^-dependence 
only if the the collinear partons are a quark and an anti-quark or two gluons with opposite 
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helicities. Therefore if quark and gluon jets could be separated experimentally (a difficult 
task) one could get a larger (cos2</?) signal by only considering collinear gluon pairs or 
quark anti-quark pairs. 

This method provides a distinctive probe of tr(G 3 ). The cos2<^ expectation value is 
dominated by the tr(G 3 ) correction and does not receive contributions from the tree-level 
five-parton QCD cross-section or from dimension-6 operators that can be related to four- 
quark contact operators. One source of background to this measurement is higher-order 
contributions to the QCD cross-section which have a non-trivial dependence on azimuthal 
rotations in the collinear region. As discussed earlier, these contributions are probably 
quite small. Some azimuthal dependence could also arise from color-strings associated 
with the hadronization of the out-going partons. The magnitude of these effects can be 
estimated by combining the exact five-parton cross-section with Monte-Carlo simulations 
for parton showers and hadronization, but we have not done so here. Also, the different de- 
pendence of the effective operator contribution on the parton center-of-mass energy could 
be used to separate its signal from these QCD backgrounds. In addition to these theo- 
retical issues, there are experimental difficulties associated with measuring the azimuthal 
dependence of almost collinear jets. In particular, both the algorithm used to define these 
jets and the detector should not induce a significant spurious azimuthal dependence; the 
cos 2<p asymmetry is only a few percent for the value of A considered here. 

So far we have considered the effects of tr(G 3 ) in two regions: one where the three 
jets are well separated, and the other where two jets are almost collinear. The tr(G 3 ) 
correction is positive in the first region and negative in the second. (The sign of the effect 
is reversed if the coefficient of the operator is negative.) One can use this fact to enhance 
the tr(G 3 ) signal by studying normalized distributions in a variable which parametrizes 
the extent of collinearity, such as the angle between jets 4 and 5 or X3, over a wide region 
which includes both the region where the jets are well separated and the region where they 
become collinear. Since the tr(G 3 ) corrections change sign over this region the resulting 
normalized distributions are more sensitive to the tr(G 3 ) signal compared to studies which 
avoid the collinear region altogether. This is indeed the case as one can see in Fig. 4 where 
we plot the normalized distributions in X3 with X3 varying from its minimum value, 2/3, 
all the way up to 0.95. Here we combine the cuts on cos6> and on ip of set A (26) which 
ensure that the five partons are well separated outside the collinear region, with the cuts 
on x$ and on the pr of jets 4 and 5 of set B (29) so that jets 4 and 5 are still distinguishable 
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when the collinear region is approached. We label these cuts as set C, 



v 7 ! > 250 GeV , 



x 3 < 0.95 , 



£5 > 0.3 , 



X4 sin 6*34 > 



5GeV 



(32) 




cos#| < 0.8 , 
30° < ip < 150° . 



The deviation due to tr(G 3 ) is more pronounced here than in Fig. 2, where the xs 



distribution is evaluated in the region where the jets are well separated. Again, the errors 
bars shown are the estimated statistical errors assuming an integrated luminosity of 15 
pb -1 . Except in the cross-over region where it changes sign, the deviation from QCD is 
significantly larger than these errors. For x 3 between 0.7 and 0.8 the deviation from QCD 
is about twenty-five percent of the QCD result. 

To exhibit the energy dependence of the tr(G 3 ) signal we also study the difference 
between the cross-section above the cross-over point, which is roughly xs = .88, and the 
cross-section below this point, with and without the tr(G 3 ) correction. This quantity is 
sensitive to the tr(C 3 ) signal because the tr(C 3 ) correction changes sign near x 3 = .88. In 
Fig. 5 we show the difference between the sum of the last three bins in the x 3 distributions 
of Fig. 4 (bins 10, 11 and 12) and the sum of the three bins below them (bins 7, 8 and 9) 
divided by the sum of all six bins for different values of y/I, 



The deviations from QCD are very large in this case, and the different energy dependence 
of the two signals is clearly seen. However, including the form-factor (27) flattens the 
tr(G 3 ) correction and reduces the deviation from QCD to the point where it is smaller 
than, or on the order of the statistical errors. 

We have discussed three approaches to probing tr(G 3 ): one which is limited to the 
region where the three jets are well separated, one which focuses on the collinear region 
and one which combines these two regions. The last two are more sensitive to the tr(G 3 ) 
signal. All the distributions we study are normalized distributions. This should reduce 
some of the systematic errors in performing the measurements. 




(33) 



23 



In most of our analysis we have included the effects of a form-factor on the deviations 
from QCD due to the operator tr(G 3 ). This is not usually done in similar searches for quark 
substructure. Studying the effects of a form-factor on the tr(G 3 ) signal became necessary 
in the context of dijet production, where this operator only contributes at order (9(1/A 4 ) so 
that the parametrization in terms of the effective operator breaks down at energies where 
the signal becomes appreciable. In that case the form-factor greatly reduces the tr(G 3 ) 
signal [6]. In the case at hand the effect of the form- factor is smaller, but non-negligible. 
Furthermore, unlike similar searches for quark substructure which use the pr or some 
other energy-dependent distributions, and are therefore sensitive to the effects of form- 
factors, most of the distributions we consider here are distributions in energy-independent 
quantities such as angles and energy fractions. A form-factor therefore only reduces the 
deviation from QCD uniformly over the entire range, by a factor which is equal to the 
energy averaged value of the form-factor, without changing the qualitative shape of the 
distributions. 

Throughout this discussion, we took the cut on the minimum center-of-mass energy 
of the jets, \/§, to be rather low - 250 GeV, as used in early CDF studies [23] of three-jet 
events. With the larger integrated luminosity now available at the Tevatron, this cut can 
be increased in order to enhance the tr(G 3 ) signal. Varying this cut also allows one to 
distinguish between the tr(G 3 ) signal and higher order QCD effects which have different 
power-law dependences on the parton center-of-mass energy 

Finally, a remark regarding the choice of of the coupling constant a s (p) is in order. 
In regions where the three jets are well defined the results are not very sensitive to this 
choice but it becomes relevant in regions where two of the jets approach collinearity. The 
cos2</? expectation value in this region is obtained by normalizing the cos2</? weighted 
distributions by the non-weighted distributions so it is insensitive to the value of the 
coupling. However, the xs distribution which interpolates between the collinear and the 
non-collinear regions would probably be affected. Since in this paper we only present 
qualitative tree- level results this effect was not taken into account. 

5. Conclusions 

As discussed in section 1, the operator tr(G 3 ) cannot be reliably probed through dijet 
production or inclusive jet production in hadron colliders. Jet production in e + e~ machines 
cannot yield good bounds on the scale A because of the low characteristic center-of-mass 
energy. This leaves three jet production in hadron colliders as the best place to probe the 
gluonic operator tr(G 3 ). In this case tr(G 3 ) contributes at leading order in 1/A 2 . 
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We have suggested two complementary ways of detecting the tr(G 3 ) signal. One 
involves studying the region where two of the jets are almost collinear, and the other 
interpolates between the collinear region and the region where the three jets are well 
separated. These approaches appear to have better analyzing power than tests that avoid 
the collinear region altogether, since the deviations from QCD are more pronounced in the 
collinear region. The two approaches may still be rather challenging experimentally. 

In the collinear region the tr(G 3 ) correction, unlike the tree- level QCD cross-section, 
has a non-trivial behavior under azimuthal rotations of the collinear momenta that leave 
their sum unchanged. Weighting the cross-section by the appropriate function of this 
azimuthal angle can in principle allow us to separate the tr(G 3 ) signal from QCD. The 
fact that the tree-level QCD cross-section is invariant under these azimuthal rotations 
follows from the helicity structure of the four-parton QCD amplitudes at tree-level, which 
in turn follows from the supersymmetric Ward identities. Thus, measuring the dependence 
of the collinear three-jet cross-section on this azimuthal angle provides a test of the helicity 
structure of tree-level QCD and the SWI, as well as a method for probing new physics in 
the strong-interaction sector via the effective operator tr(G 3 ). 
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Appendix: Collinear properties of color-ordered sub-amplitudes 

In this appendix we discuss additional details of color ordered sub-amplitudes, the 
helicity basis and collinear factorization. Except for the behavior of the tr(G 3 ) three-gluon 
vertex in the collinear limits, all the results mentioned are well known (see the review 
paper by Mangano and Parke [16] and references therein) and are summarized here for the 
reader's convenience. 

The gluon sub-amplitudes (2) have the following properties: 

• A n (l, 2, . . . , n) is invariant under cyclic permutations. 

• A n {n, n - 1, . . . , 1) = (-l) n A n (l, 2, . . . , n) 

• A n (l, 2, 3, . . . , n) + A n {2, 1, 3, . . . , n) + ■ ■ ■ + A n (2, 3, . . . , 1, n) = 
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The last equation is due to the decoupling of the U(l) gauge boson in a U(N) gauge 
theory. These relations reduce the number of independent sub-amplitudes that have to be 
calculated. 

We now turn to discuss collinear factorization in the color-ordered basis. Equa- 
tion (6) of section 2 can be proven from the string representation of gauge theory tree 
amplitudes [16] but here we will briefly discuss it from the point of view of Feynman di- 
agrams. In the limit when partons % and j become collinear, the n-parton color ordered 
sub-amplitude with n > 5 can become singular as the relevant propagator goes on-shell. 
This only happens if the collinear partons are color- adjacent; they have to be attached 
to the same three-parton vertex to give rise to an on-shell propagator. The singularity 
may be removed by powers of the Lorentz invariant s^- coming from the vertex. However, 
for most helicity choices of i and j and A, the QCD three-parton vertex contracted with 
the polarization vectors of i and j (or spinors, for quarks), behaves as y/Si] at leading 
order so that the singularity remains, although it is reduced to 1/^/SiJ. Furthermore, the 
coefficient of the ^phj term in the vertex can be replaced, up to a numerical factor, by the 
polarization vector of an effective external parton with momentum p = ki + kj and helicity 
A. The remaining part of the diagram, contracted with the effective parton, becomes an 
(n — 1) -parton diagram, and summing over all possible diagrams gives the (n — 1) -parton 
amplitude, multiplied by the splitting function. 

We now show that no additional collinear singularities arise from the gluon operator 
tr(G 3 ). The only issue is for the quark-free g — > gg singularity. Contracting the new three 
gluon vertex with the polarization vectors e\ and e 2 of the external gluons 1 and 2 and 
using momentum conservation and the gluons' transversality, e\ • k\ = s 2 • k 2 = 0, this 
vertex can be simplified to the form 



up to an overall factor which does not depend on the momenta fci, k 2 of the gluonsf. 

In the limit where k\ and k<i become collinear, we can write k\ = zp + q, k 2 = 
(1 — z)p — q, where z is a number and g M — > 0. Eliminating p from these relations we can 
express k\ as a combination of q and k 2 and vice versa. Using e\ • k\ = £2 • ^2 = the first 
term in eqn. (34) becomes 



j If the two gluons have opposite helicities, one can choose the reference momentum of 1 to be &2 and 
vice versa so that the vertex vanishes for arbitrary fciand ki. 




(34) 



ki ■ e 2 k 2 ■ ei oc q ■ e 2 q ■ £1 • 
But since the on-shell propagator is l/(/ci + k 2 ) 2 and 

(h + h) 2 = 2h-k 2 = p 2 = 0(q 2 ) = O(p-q) , 
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we see that both terms in the vertex go to zero at least as fast as {k\ + k 2 ) 2 so that the 
singularity cancels. 

The explicit representation of the spinor products that we use follows ref. [16]. For 
two momenta with positive energies: 

(ij) = y/s~ expiicpij) , (35) 

where = 2pi ■ pj , and 

cos fa = (PiPj -PjPt)/^ptppij , sin (Pi j = (p 2 p+ - p 2 p+ ) / pfp+ Slj , 

where p ± = (p° ±p 3 ). If any of the momenta has negative energy, the spinor product is 
calculated with minus that momentum and then multiplied by i for each negative energy 
momentum, [ij] can be calculated using (ij) [j i] = Sij and the antisymmetry of the spinor 
products, (ij) = - (ji), [ij] = - [ji]. 

We collect here the tree- level splitting- functions used to infer the tr(G 3 ) correction to 
five-parton amplitudes. 

The g — > gg splitting functions are 



Split + 
Split. 

Split. 
Split + 
Split. 
Split + 
Split + 
Split. 



* ,3 



» ,3 



« ,3 



+ 



i ,3 



0, 
0, 





z) (ij)' 




1 




-z) [ij] 
z 2 




~ z) [ij] 


(1- 


zf 


z 


z) (ij)' 

2 


y/z(l- 


z) {ij)' 


(1 


-zf 




~z) [ij] 



(36) 
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the g — > qq splitting functions are 

Split + (?+ <T) 

Split + (q-,q+) 
Sp\it_(q+,q-) 
Split. (q~,q + ) 



z l/2 (1 _^)3/2 

yfzQ - z) (qq} 1 
z 3/2( 1 _ z y/2 

y/z(l-z) (qq) 

z 3/2 (1 _ z) l/2 

^/zJT^z) [q q] ' 
z i/2 (1 _^)3/2 

^/z(T^~z) [q q] ' 
and the q ^ qg and q ^ gq splitting functions are 

^3/2 

Split + (g 

Spiit + (g",r) 

Split. 

Split_ , i~) 
Split_(z+, g+) 
Split _ (z~, g + ) 
Split + (i + ,g~) 

Spiit + (r,Q _ ) 



y/z(l-z) (qi)' 
z 1 ' 2 
y/z(l-z) [qi]' 
z 1 / 2 

y/z(l-z) {qi}' 
z 3 ' 2 
y/z(l-z) [qi]' 
(1-z) 1 / 2 
yfzjl - z) (iq) 7 
(1-z) 3 / 2 
\/z{l - z) [i q] 7 
(1-z) 3 / 2 
y/z(l - z) (iq) 7 
(1-z) 1 / 2 
y/z(l - z) [iq] 



(37) 



(38) 



Here z denotes the collinear momentum fraction, so that ki = zp, kj = (1 — z)p, with p 
denoting the sum of the collinear momenta. Note that reversing all helicities (parity) is 
equivalent to exchanging ( ) and [ ] , and inserting a minus sign. 
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Figure Captions: 



Figure 1: The three jet event in the parton center-of-mass frame. The parton momenta 
are not drawn to scale. 

Figure 2: Effect of tr(G 3 ) on various distributions in three jet event variables at ^/s = 
1.8 TeV in the region where the jets are well separated. The solid line gives the QCD 
prediction, and the dashed line includes the tr(G 3 ) correction with A = 1 TeV, with the 
cuts of set A (equation (26)). The error bars are the estimated statistical errors for an 
integrated luminosity of 15 pb -1 . 

Figure 3: Effect of tr(G 3 ) on various distributions of three jet event variables at ^/s = 
1.8 TeV in the region where two of the jets are almost collinear. The solid line gives the 
QCD prediction, the other two lines include the tr(G 3 ) correction with A = 1 TeV without 
(dashes) and with (dot-dashes) the form-factor (1 + 47rs/(3A 2 )) 2 . Here we impose the 
cuts of set B (equation (29)). The error bars are the estimated statistical errors for an 
integrated luminosity of 15 pb -1 . 

Figure 4: Effect of tr(G 3 ) on the x s distribution at yfs = 1.8 TeV in the cross-over region. 
The solid line gives the QCD prediction, the other two lines include the tr(G 3 ) correction 
with A = 1 TeV without (dashes) and with (dot-dashes) the form-factor (1 + 4rrs/ (3A 2 )) 2 . 
Here we impose the cuts of set C (equation (32)). The error bars are the estimated 
statistical errors for an integrated luminosity of 15 pb -1 . 

Figure 5: The difference between the last three (10,11,12) bins and the three bins below 
them (7,8,9) of the X3 distribution of Fig. 3, divided by the sum of the bins, as a function 
of v 7 ^. The error bars are the estimated statistical errors for an integrated luminosity of 
15 pb" 1 . 
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